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SUMMARY 


Two theoretical methods are presented for optimizing multielement air- 
foils to obtain maximum lift. The analyses assume that the shapes of the 
various high-lift elements are fixed. The objective of the design procedures 
is then to determine the optimum location and/or deflection of the leading- 
and trailing-edge devices. The first analysis determines the optimum hori- 
zontal and vertical location and the deflection of a leading-edge slat. The 
structure of the flow field is calculated by iteratively coupling potential- 
flow and boundary-layer analysis. This design procedure does not require 
that flow separation effects be modeled. The second analysis determines the 
slat and flap deflection required to maximize the lift of a three-element 
airfoil. This approach requires that the effects of flow separation from 
one or more of the airfoil elements be taken into account. The theoretical 
results are in good agreement with results of a wind-tunnel test used to 
corroborate the predicted optimum slat and flap positions. 


INTRODUCTION 


To achieve acceptable takeoff and landing performance, modem fixed- 
wing aircraft typically use mechanical high-lift devices such as leading-edge 
slats and slotted flaps. An essential step in the development of high-lift 
wings is the determination of the position (horizontal location, vertical 
location, and deflection) of the slat and/or flaps required to maximize 
aerodynamic performance. Theoretical methods which predict the optimum 
position of slats and flaps for multicomponent airfoils are being developed 
to; (1) reduce the number of high- lift configurations which must be evaluated 
experimentally; (2) reduce substantially the amount of wind-tunnel testing 
required to optimize a multicomponent airfoil; and (3) Insure that experi- 
mentally determined optimums are true optimums. Two criteria generally used 
to evaluate the performance of a high-lift airfoil are maximum lift coef- 
ficient and maximum lift to drag ratio (Cj^/Cd)^^^ a given lift 

coefficient. This paper addresses the problem of optimization for maximum 
lift. At present, Jrag is not considered a criterion in the theoretical 
optimization method because the drag of multielement airfoils cannot be 
predicted accurately. 
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Configuration optimization of multiele<aent airfoils for lift 

is separable into two general classes of problems. This classification 
depends on the level of sophistication required of the aerodynamic analysis 
to determine the structure of the flow field and to evaluate the airfoil 
performance. Configurations which can be optimized without modeling the 

^ problems. Configuration 

optimization which requires that the effects of flow separation be modeled 
are called Class 2 problems. For Class 1 optimization, viscous effects are 
determined by an analysis based on laminar and turbulent boundary-layer 
flows on the various airfoil surfaces. This type of analysis is simpler, 
requires less computation effort, and is generally more reliable than the 
separated flow modeling techniques required for Class 2 problems. 

Section I outlines the aerodynamic analysis including a brief description 
of the separated flow model. The Class 1 design procedure is then used in 
action II to determine the optimum position of a leading-edge slat. 
Optimization of the deflection angles of a slat and of a slotted flap for 
a three-element configuration is presented as a Class 2 design problem in 
Section III. Section IV briefly describes the experimental study conducted 
to obtain data for verification of the analysis techniques presented in this 

paper, and Section V compares theo y with experiment for both optimization 
niGtiinods • 


SYMBOLS 

A wing area ■ cb 

fR wing aspect ratio - b^/A 

b wing span 

c airfoil reference chord » 1.7 m 



C 

P 

P 

Re 


section drag coefficient ■ 

(1/2) PiUoo^c 

skin-friction coefficient at the upper surface trailing edge 
of airfoil component 

airfoil lift coefficient - 

(1/2) p U 2c 
00 00 

wing lift coefficient = 

(1/2) P„U^2A 

p-p 

airfoil surface pressure coefficient - t—. 

(1/2) pU 2 
00 

static pressure 

U c 

reference chord Reynolds number ■ 


« 
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00 
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x,y 

a 

6 

V 
P 
T 


airfoil surface (measured from stagnation 

dimensionless distance from separation point (eq. (6)) 
s-dlrectlon velocity component at out?r edge of boundary layer 
free-stream velocity 

velocity component normal to airfoil surface 

horizontal and vertical coordinates in coordinate system of main 
airfoil 

angle of attack, deg 
component deflection, deg 
kinematic viscosity 
fluid density 
wall shear stress 


Subscripts : 


m 


flap 

main element 
slat 

free stream 


I. AERODYNAMIC ANALYSIS 


General Comments 

of . a»1yl8 18 detormln. th. po.ltloo 

tn element (or elements) required to achieve maximum lift. Up 

f” '■'’“‘“'I 1“ ’P^lfy position of o.ch elemeoi. 
^686 po8ltion variables plus the angle of attack of the overall conflgura- 

f ®“lti‘*i®enslon design space which is searched to determine 
the optlmm position. As a result, not only must the ae -odynamic analysis 

tlonr^hif flow-field structure for a wide Variety of conflgura- 

tlons, but it must also be computationally efficient if the optimization 
technique is to be a useful design tool. It should be noted that the 
prlM^ objective of the theoretical analysis is to determine the position 
high-lift element (or elements) required to achieve maximum lift. 


Although the present analysis should give a reasonable estimate of the 
maximum lift coefficient of an optimized configuration, a more accurate 
computation of will be the subject of future study. 


Potential-Flow Calculation 

The potent ial-f low calculation is a singularity method where each of 
the airfoil contours (fig. 1) is represented by a large number (from 40 to 
80) of straight line segments. A linear distribution of vorticity is placed 
on each of these segments. The present potential-flow analysis (ref. 1) 
not only provides the accurate surface pressure distributions required for 
subsequent boundary-layer analyses, but it is also computationally efficient. 
This efficiency is essential because a typical optimization calculation 
requires from 200 to 500 potential-flow solutions to complete the design 
procedure. 


Bound ary- Layer Analysis 

The structure of the laminar, transitional, and turbulent boundary layers 
is calculated using the pressure field determined from the potential-flow 
analysis. For the three-element configuration considered in this study, a 
complete optimization calculation requires between 500 and 2000 solutions 
for conventional boundary layers such as those on the upper and lower surface 
of the leading-edge slat shown in figure 1. The need for computational 
efficiency resulted in the selection of the finite-difference method of 
Blottner (ref. 2) for solving the boundary-layer equations. The correlations 
of Smith (ref. 3) are used to determine the point of laminar instability. 

The transitional and turbulent b.^uni ary- layer eddy viscosity model of 
Cebeci (ref. 4) and Cebeci, Kaups, Mosinakis and Rehn (ref. 5) is then used 
to calculate Reynolds stresses. If laminar boundary-layer separation occurs, 
the empirical criteria of Caster (ref. 6) is used to determine whether 
turbulent reattachment occurs. If reattachment is indicated, the calculation 
is continued as a turbulent boundary layer. If the criterion of Caster 
indicates catastrophic flow separation, the boundary-layer calculation is 
terminated. A boundary-layer calculation of this type takes only 0.1 to 
0.2 sec on the CDC 7600 computer. 

The structure of the confluent boundary layer on the flap upper surface 
(fig. 1) is determined using the finite difference method of referenc* 7, 
a method which takes from 20 to 30 sec on the CDC 7600 computer for e. ^h 
confluent-boundary- layer calculation. 


Potential-Flow/Boundary-Layer Coupling 

The effect of boundary-layer displacement and entrainment on the inviscid 

flow is simulated by linear source distributions (fig. 1) placed on the 

straight line segments used to define the airfoil contours. The strength 

of the source distribution is equal to d(u 6*)/ds where u is the velocity 

e e 


at the edge of the boundary layer, 6* Is the local boundary-layer dlsplace- 
^nt thickness, and s is the arc length measured along the airfoil surface. 
This type of potential-flow/boundary-layer coupling is particularly useful 
when optimizing the position of high-lift devices because the displacement 
effect of any particular boundary layer can conveniently be "frozen." For 
example, the boundary layer on the lower surface of the flap (fig. 1) is 
Mt appreciably affected by changes in the leading-edge slat position. 
Therefore, the flap lower-surface boundary layer need only be computed for 
the initial slat position. The displacement effect for subsequent slat 
positions is then based on the flap lower-surface source strength distribution 
as calculated with the slate in the initial position. 


Separated Flow Analysis 

The separated flow model developed for the present study is similar 
to the source distribution method of Jacob and Steinbach (ref. 8). The 
separated flow model requires that the static pressure distributions in 
the separated zone (fig. 2) be nearly constant. This constant pressure 
condition in the separated region is approximated by requiring that the 
static pressure at (the central panel in the separated zone) and the 
the static pressure at Sj (the trailing edge of the airfoil) both be equal 
the static pressure at s^ (the boundary- layer separation location). 


C (s ) - C (s ) (1) 

P 2 Pi ' ^ 

and 


C (8 ) - C (a ) 
P 3 Pi 


( 2 ) 


A two-parameter source distribution located on the surface of the airfoil 
in the separated zone (fig. 3) is introduced to enforce the above conditions. 
The functional form of the source distribution Is given by 


V » 2.5 V s* 
max 


0- S*- 0.4 


(3) 


V - V 


max 


0.4-^,. -0.6 


(4) 


V - (v - V )(5s*-3)2/4 
max max 3 


0.6- s*- 1.0 


(5) 


where 


8* - (s - Sj)/(S 3 - Sj). 


( 6 ) 


The separation location calculated from a previous iteration is used to 

determine the maximum source strength in the separated zone (v ) and the 

max 
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source strength at the airfoil trailing edge (v^) so that equations (1) 
and (2) are satisfied. The airfoil static pressure distribution obtained 
from this solution is then used to recompute the airfoil upper-surface 
boundary layer and a new separation point is determined. This process is 
repeated until convergence is achieved. 

Figure 4 coiiq>ares the theoretically predicted lift for a range of 
angle of attack with experiment (ref. 9) for a GA(W)-1 airfoil. Figure 5 
presents a similar comparison for a three-element configuration tested b; 
Foster, Irwin, and Williams (ref. 10). The agreement between theory and 
experiment is good for both cas 2 s at angles of attack below stall. At angles 
of attack near stall the lift is underpredicted by approximately 8% for the 
GA(W)-1 airfoil and by 3% for the three-element configuration. 


II. SLAT OPTIMIZATION: A CLASS 1 PROBLEM 


The objective of this optimization procedure is to theoretically 
determine the position (horizontal location, vertical location, ard deflec- 
tion) of a leading-edge slat for maximum lift based on aerodynamic calcula- 
tions which do not model the effect of flow separation. The analysis is 
based on the premise that, at the maximum lift coefficient at which 
attached flow can be maintained on the main airfoil, the optimum slat 
position for maximum lift minimizes the suction peak on the upper surface 
of the main element. For configurations which stall abruptly this is a 
reasonable assumption. However, for configurations which stall gradually, 
the validity of such an assumption remains to be proven. 

The constrained function minimization method described by Vanderplaats 
in reference 11 and Vanderplaats and Moses (ref. 12) is used to numerically 
optimize the slat position so that the suction peak in the leading-edge 
region on the main component is minimized. This minimization is performed 
with the main element at a fixed angle of attack. The design variables 
(fig. 6) are slat horizontal and vertical location, and slat deflection. 

All three parameters are referenced to the slat trailing edge and can be 
varied independently. 

Slat translation and deflection are subject to two constraints. The 
first constraint - that there be no flow separation on the slat upper 
surface -is applied to prevent the slat from moving into a separated flow 
regime. Numerically, this condition is satisfied by requiring that the 
skin-friction coefficient at the trailing edge of the upper surface of the 
slate, Cf , be slightly positive (i.e., Cf i 0.0001). The second con- 
straint riquires that the slat trailing edge be no closer than 0.02 c to 
the surface of main element (a reasonable slat gap for this configuration 
(ref. 10)). This constraint prevents the slat wake from merging with the 
boundary layer on the wing upper surface so that the computationally time 
consuming confluent boundary-layer calculation for the viscous flow above 
the main element can be replaced by a computationally efficient conventional 
boundary- layer calculation. 
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III. SLAT/FLAP DEFLECTION OPTIMIZATION: 
A CLASS 2 PROBLEM 
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between the wakes and the upper-surface boundary layers. It was also 
assumed, therefore, that the viscous flows on the upper surface of the main 
wing and flap could be analyzed as conventional boundary layers. 


IV. EXPERIMENTAL OPTIMIZATION 


As part of this study tests were conducted In the NASA Ames 40-by-80- 
Foot Wind Tunnel to experimentally determine optimum slat and flap positions 
for comparison with the theoretical predictions. The rectangular planform 
wing used In these tests (fig. 13) Is equipped with a full-span leading-edge 
slat and a full-span single-slotted flap. The slat chord Is 0.17 c and the 
flap chord Is 0.40 c. The basic airfoil section Is on RAE 2815. A detailed 
description of the slat, main airfoil, and flap shapes can be found In 
reference 10. The wing span Is 16 m and the extended chord Is 2.15 m. The 
relatively high aspect ratio of 7.5 and the rectangular planform result in 
a configuration with nearly two-dimensional flow over much of the wing span. 
Also, the use of a high aspect-ratio finite wing eliminates adverse wind- 
tunnel wall interference effects associated with two-dimensional high-lift 
airfoils which span the entire test section. 

The slat and flap brackets permit continuous adjustment of the hori- 
zontal and vertical locations and the deflection of both elements. Data 
taken Included lift using the wind-tunnel balance and surface static 
pressures along the model centerline. All data were obtained at a Reynolds 
number of 3.8x10® and a Mach number of 0.10. 


V. COMPARISON OF THEORY AND EXPERIMENT 


The experimentally and theoretically optimized slat positions are com- 
pared in figure 14 for a fixed flap deflection of 10°. The theoretical 
prediction Is based on the Class 1 design procedure presented In Section II. 

It Ig to be noted that the largest differences uatween theory and experiment 
Is 4 of slat deflection. The measured wing maximum lift coefficient with 
the slat In the theoretically optimized position differs only 4% from the 
experimentally measured maximum lift coefficient with the slat In the experi- 
mentally optimized position. 

The comparison between theory and experiment for the slat/flap deflection 
optimization is summarized In figures 11 and 12. The experimentally measured 
maximum lift coefficient of the wing, CLmax> presented together with the 
wing centerline section maximum lift coefficient, Ci . The experimental 
value of C£ was determined by correcting the meSlilred C^ for aspect 
ratio effects using the lifting line theory described by Glauert^(ref . 13). 

The ekperimentally determined slat and flap deflections required to maxi- 
mir.'' the total wing lift are In good agreement with the optimum slat and 
flap deflections determined from the Class 2 analysis presented In Section III. 
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VI. CONCLUDING REMARKS 
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CONVENTIONAL CONFLUENT CONFLJENT 

BOUNDARY LAYER BOUNDARY LAYER BOUNDARY LAYER 



PIECEWISE LINEAR VORTICITY AND 
SOURCES DISTRIBUTED ON AIRFOIL 
SURFACES 


Figure 1.- Flow about multielement airfoil. 



*2 " CENTER PANEL IN SEPARATED ZONE 
I3 » TRAILING EDGE OF AIRFOIL 

REQUIRE: Cp(it ) - Cp(at ij) - Cp(«t •,) 


Figure 2.- Separated flow model 



SOURCE STRENGTH. 



Figure 3.- Two-parameter source distribution in separated zone. 



VISCOUS THEORY 

INVISCIO THEORY 

EXPERIMENTAL 
(REr. 9) 


Figure 4.- Comparison of experimental and calculated lift on a GA(W)-1 
airfoil. Re - 5.7 x in® : M - 0.15. 
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Figure 


Figure 6,- 



INVISCIO THEORY 
— VISCOUS THEORY 
O EXPERIMENT 


5.- Comparison of experimental and calculated lift on a three- 
element airfoil. Re - 3.8 X io6; M - 0.20. 
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X, = x-COORDINATE OF SLAT TRAILING EDGE 
V, • V-COOROINATE OF SLAT TRAILING EDGE 
S, • SLAT DEFLECTION / VARIABLES 


Aerodynamic analysis for leading-edge 
for maximum lift. 


slat-posltlon optimization 


iTERATIOM 



Figure 7.- Minimization of main airfoil suction pea«. at an angle of attack 
of 15 . « 10®; Re • 3.8 x lo*; M « 0.10. 



Figure 8.- Pressure distribution on a three-element airfoil with ?lat 
positioned to minimize the suction peak on the main airfoil. 

6f - 10®; Re - 3.8 x lo*: M - 0.10. 
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ANGLE OF ATTACK, a, deg 

Figure 9.- Effect of configuration angle of attack on the skin-friction 
coefficient at the trailing edge of the main component and on the slat 
position required to minimize the suction peak on the main component. 

= 10®; Re - 3.8 x 10®; M = 0.10. 
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BOUNDARY LAYER BOUNDARY LAYER 


CONVENTIONAL 
BOUNDARY LAYER 



ft, ■ SLAT DEFLECTION 1 DESIGN 
ft, > FLAP DEFLECTION I VARIABLES 


Figure 10.- Aerodynamic analysis for leading-edge slat and trailing-edge 
flap-deflection optimization for maximum lift. 
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SLAT DEFLECTION, dig 

deflection on the computed 

maximum airfoil lift with the effect of slat deflection on th( 

measured wing maximum lift" Ct for three-element configuration. 

Re = 3.8 X 10®; M = 0.10. 
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FLAP DEFLECTION. S|. dig 


Figure 12.- Comparison of the effect of flao deflecUnn «« .-k.. 
R. . M ‘ '""'-•‘P""'' configuration. 



Figure 13*- Wind-tunnel model used to determine experimentally optimum 
slat and flap positions for maximum lift. 


SLAT LOCATION 

x,/c 



(EXPERIMENTAL) 

A 

-0.007 

0.017 

-36 

3.16 

B 

-0.009 

0.010 

-40 

3.28 



Figure 14.- Comparison of theoretically and experimentally optimized 
leading-edge slat position. Position A is theoretical; position 
B is experimental; 6^ = 10^; Re =* 3.8 x 10^. 




